meson production in Au + Au and p + p collisions at ^/s^ = 200 GeV 
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We report the STAR measurement of (j> meson production in Au + Au and p + p collisions at 
y'Sj^j^ = 200 GeV. Using the event mixing technique, the cp spectra and yields are obtained at mid- 
rapidity for five centrality bins in Au+Au collisions and for non-singly-diffractive p+p collisions. It is 
found that the (j) transverse momentum distributions from Au + Au collisions are better fitted with a 
single-exponential while the p-l-p spectrum is better described by a double-exponential distribution. 
The measured nuclear modification factors indicate that production in central Au + Au collisions is 
suppressed relative to peripheral collisions when scaled by the number of binary collisions {{Ntin})- 
The systematics of (pt) versus centrality and the constant 4>/K~ ratio versus beam species, centrality, 
and collision energy rule out kaon coalescence as the dominant mechanism for (j) production. 

PACS numbers: 25.75.Dw 

In elementary collisions the production of the (j) me- pectations from p + p collisions 0, la, ISl- Theoretical 
son, the lightest bound state of strange quarks {ss) , is calculations have tried to address the origins of this en- 
suppressed because of the OZI rule 0, SH- ^"^ heavy- hancement 0, S 0I- The 4> meson is also thought to 
ion collisions, however, strange quarks are produced co- have a small hadronic cross-section |lO| and may provide 
piously and 4> enhancement is observed relative to ex- direct information about the dense matter at hadron for- 



mation without perturbations from co-moving hadrons. 
For these reasons, </> production in relativistic nuclear col- 
hsions has been of great interest. 

The mechanism for (j) production in high energy col- 
hsions has remained an open issue. A naive interpreta- 
tion of the (j) enhancement observed in heavy-ion coUi- 
sions would be that the (f) is produced hadronically via 
KK — > (f). Hadronic rescattering models such as RQMD 
and UrQMD [lll[l2|j implementing such processes, pre- 
dict an increase in the (t>/K~ ratio as a function of the 
number of participants. Rescattering models also predict 
similar increases in the ipt) of the proton and <j) meson. 

The nuclear modification factors {Raa and Rcp) of 
the (f> meson are important in differentiating between 
mass and particle species ordering. Current measure- 
ments of identified hadrons by STAR (A and K 'a) s how 
that Rcp for the A differs from that of the if gill. It 
is difficult, however, to determine whether this difference 
is related to the mass of the particle or the type of the 
particle (whether it is a baryon or a meson) since there is 
a significant mass difference between the A and the Kg. 
The (p, however, has a mass that is similar to that of the 
A, yet is a meson. A direct comparison of the (p Rcp and 
Raa with these previous measurements will provide more 
insight into this mass vs. particle species dependence. 

The STAR detector |lj| consists of several sub-systems 
in a large solenoidal analyzing magnet. For the data 
taken during the second RHfC run (2001-2002) presented 
here, the experimental setup consisted of a Time Projec- 
tion Chamber (TPC), a Central Trigger Barrel (CTB), a 
pair of Beam-Beam Counters (BBC), and two Zero De- 
gree Calorimeters (ZDC). The ZDC's are used as the ex- 
perimental trigger for Au + Au collisions while the BBC's 
are used for the p + p trigger. 

The results presented here were obtained from about 
2.1 million minimum-bias Au+Au events, 1.1 million cen- 
tral Au-\- Au events and 6.5 million non-singly-diffractive 
(NSD) p+p events. Reconstruction of the </> was accom- 
plished by calculating the invariant mass {rriinv), trans- 
verse momentum (pt), and rapidity (y) of pairs that 
formed from all permutations of candidate K^ with K~. 
The resulting rriiny distribution consisted of the (p signal 
atop a large background that is predominantly combina- 
torial. The shape of the combinatorial background was 
calculated using the mixed-event technique [T^ [ia| ■ 

For the centrality measurement, the raw hadron mul- 
tiplicity distribution within a pseudo-rapidity window 
1 77 1 < 0.5 is divided into five bins corresponding to 50- 
80%, 30-50%, 10-30%, top 10% and top 5% of the mea- 
sured cross section for Au + Au collisions. Events are 
selected with a primary vertex z position from the center 
of the TPC of |z| < 25 cm for Au + Au collisions and 
\z\ < 50 cm for p + p collisions, where z is along the 
beam axis. These events are further divided according to 
z to reduce acceptance-induced distortions in the mixed- 
event background. Correlations in the background due 
to elliptic flow were minimized by mixing events with 
similar reaction plane angles. Consistent results are ob- 



tained when we construct the background distribution 
using like-sign pairs from the same-event. 

Particle identification (PID) is achieved by correlating 
the ionization energy loss (dE/dx) of charged particles 
in the TPC gas with their measured momentum. The 
measured (dE/dx) is reasonably well described by the 
Bethe-Bloch function 10, 17] smeared with a resolution 
of width a. Tracks within 2cr of the kaon Bethe-Bloch 
curve are selected for this analysis. 
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FIG. 1: The transverse mass distributions from Au + Au 
(circles) and p + p (squares) collisions at 200 GeV. For clarity, 
some Au + Au distributions for different centralities are scaled 
by factors. The top 5% data are obtained from the central 
trigger data set. All other distributions are obtained from the 
minimum- bias data set. Dashed lines represent the exponen- 
tial fits to the distributions and the dotted-dashed line is the 
result of a double-exponential fit to the distribution from p-l-p 
collisions. Error bars are statistical errors only. 

To obtain the (p spectra, same-event and mixed-event 
distributions are accumulated and background subtrac- 
tion is done in each pt, y and centrality bin. The mixed- 
event background m.inv distribution is normalized to the 
same-event rriinv distribution in the region above the 
(p mass (1.04 < rriinv < 1-2 GeV/c^). A small, smooth 
residual background can remain near the cp peak in the 
subtracted minv distribution, because the mixed-event 
sample does not perfectly account for the production of 
background pairs (kaons and/or pions from PID leak- 
through) that are correlated, either by Coulomb or other 
interactions or by such instrumental effects as track merg- 
ing p^ . The raw yield in each bin is then determined 
by fitting the background subtracted minv distribution 
to a Breit-Wigner function plus a linear background in 



a limited mass range. The measured mass and width of 
the <j) are consistent with the value listed by the Particle 
Data Group [l3| convoluted with detector resolution. 

Using GEANT and detector response simulations, the 
data are corrected for acceptance, kaon decay and track- 
ing efficiencies to obtain the final distributions presented 
here. Figure Q shows the transverse mass distributions 
from Au + Au (circles) and NSD p+p (squares) collisions 
at 200 GeV. The spectra are obtained from the rapidity 
range jy^j < 0.5. For clarity, some Au-'r Au distributions 
for different centralities are scaled by factors indicated 
in the figure. Dashed lines represent exponential fits to 
the distributions and the dotted-dashed line represents a 
double-exponential fit to the p + p result. 

Statistical uncertainties are shown in the figure and 
the results of the fits are listed in Tabled The main con- 
tributions to the systematic uncertainty come from fit- 
ting to the K^K^ invariant-mass distribution, tracking 
and the FID efficiency calculation. Different background 
functions and normalization factors for the mixed-event 
background were used to determine the uncertainty in 
the fitting to the invariant-mass distribution and is esti- 
mated to be about 5%. The uncertainty from tracking 
and FID efficiency is estimated, by varying the tracking 
and FID cuts on the daughter tracks, to be 8%. The 
overall systematic uncertainty in the yield, dN/dy and 
{pt) is estimated to be 11%, and includes an additional 
contribution from fitting the transverse momentum dis- 
tributions. For Au + Au collisions, the inverse slope pa- 
rameters and yields are extracted from a single exponen- 
tial function fit. For p + p collisions, however, there is an 
additional component beyond a single exponential, see 
dashed-line in figure ^ The power-law shape provides 
a better fit at the higher pt region but failed at low pt . 
Double-exponential function provided a better fit so it 
was used to extract the values of dN/dy and (pt) for the 
p + p collisions. For the heavy ion resu lts, a Boltzmann 
distribution and a thermal-|-flow model [13 are also used 
to fit the data as a check of the systematic uncertainty in 
the extrapolated yield and (pt)- The systematic uncer- 
tainty is ~ 15% in the overall normalization and < 5% 
in mean pt for the p + p data, including uncertainties in 
the vertex efficiency for very low multiplicity events. 

The system-size and beam-energy dependence of (pt), 
(f>/K~ and (jy/hr are shown in figure |21 For compari- 
son, the {pt) of the p, K~ and tt" are also shown 20]. 
In the 4'lh~ ratio at y/s^ = 200 GeV shows no sig- 
nificant dependence on centrality for Au + Au collisions, 
but decreases by about 30% for p + p collisions, see open 
symbols in plot (b). As a function of energy, see plots 
(c) and open circles in plot (d), both values of {pt) and 
(p/h" ratio increase. This indicates that the production 
of 4> mesons is sensitive to the initial conditions of the 
collision. 

The general trend for p, K~ and tt" is an increase 
in {pt) as a function of centrality, which is indicative 
of an increased transverse radial flow velocity compo- 
nent to these particles' momentum distributions. The 
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0.97 ±0.02 7.70 ±0.30 

0.95 ±0.03 6.65 ±0.35 

0.97 ±0.02 3.82 ±0.19 

1.02 ±0.03 1.72 ±0.06 

0.94 ±0.02 0.48 ±0.02 

0.82 ±0.03 0.018 ±0.001 



TABLE I: Results of (f) meson inverse slope parameter, (pt), 
and dN/dy from NSD p + p and Au + Au collisions at RHIC. 
An exponential is used for the Au + Au data while a double- 
exponential fit is used for the p + p data. All values are for 
|y| < 0.5 and only statistical errors are quoted. 



(j) (pt), however, shows no signiflcant centrality depen- 
dence. This indicates that the 4> does not participate in 
the transverse radial flow as does the p, K~ and 7r~ . This 
is expected if the 4> decouples early on in the collision be- 
fore transverse radial flow is completely built up. If the 
(j) hadronic scattering cross section is much smaller than 
that of other particles, one would not expect the (j) (pt) 
distribution to be appreciably affected by any final state 
hadronic rescatterings. In contrast to these observations, 
the RQMD predictions of (pt) for kaon, proton and cj) all 
increase as functions of centrality |llL l2ll| . 

The yield ratio (j)/K^ from this analysis is constant as 
a function of centrality and species (j)+p or Au + Au). In 
fact, for collisions above the threshold for (j) production, 
the 4>/K^ ratio is essentially independent of system size, 
e+e^ to nucleus-nucleus, and energy from a few GeV up 
to 200 GeV (Figure|21(d)) ^4,5, &,J7, 22, 23, 24]. This is 
remarkable, considering that the initial conditions of an 
e+e" collision are so drastically different from Au + Au 
collisions. This observation may indicate that the ratio 
is dominated by the hadronization process. 

Rescattering models (RQMD [III, UrQMD O) pre- 
dict that about 2/3 of (j) mesons come from kaon coales- 
cence in the final state. The centrality dependence of the 
(f)/K~ ratio alone provides a serious test of the current 
rescattering models. In these models, such as UrQMD, 
rescattering channels for (p production includes KK and 
-R'-Hyperon modes and predicts an increasing 4>/ K~' ra- 
tio vs. centrality. These models also predict an increase 
in {pt) for the proton, kaon, and (f) of 40 to 50% from pe- 
ripheral to central collisions. A comparison of the data 
to these models does not support the kaon coalescence 
production mechanism for (j) mesons. 

The particle-type dependence of the nuclear modifica- 
tion factors Raa and Rcp [i3j 119 should be sensitive to 
the production dynamics and the hadronization process 
26, 27, 28, 29, 30]. Raa is the ratio of the differen- 
tial yield in a centrality class of Au + Au collisions to 
the inelastic differential cross-section m p + p collisions, 
scaled by the overlap integral Taa — {Nun) / o'inei from 
a Glauber calculation|31j. The Glauber calculation was 
performed with Gmei = 42 ± 1 mb. The inelastic differ- 
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FIG. 2: (a) (pt) vs. measured number of charged hadrons 
(Nch) within I77I < 0.5 at 200 GeV. For comparison, the val- 
ues of (pt) for negative pions, kaons, and anti-protons are 
also shown; (b) Ratios of N{(j))/'N(K^), filled symbols, and 
N(^)/N(/i~), open symbols, vs. Nch; (c) (pt) vs. center- 
of-mass beam energy from central nucleus-nucleus (filled cir- 
cles) and p + p collisions (filled triangles); (d) Ratios of 
N ((/!>) /N(i('~) from central nucleus- nucleus collisions, filled 
circles, and N((^)/N(/i~), open circles, vs. center-of-mass 
beam energy. N(<^)/N(_R'~) ratio from e'^e~ collisions (open 
squares) are also shown. Note: All plots are from mid- 
rapidity. Both the statistical and systematic errors are shown 
for the 200 GeV STAR data, while only statistical errors are 
shown for the energy dependence of the particle ratios. 
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FIG. 3: Rcp (a): The ratio of central (top 5%) over pe- 
ripheral (60-80%) (Rcp) normalized by (Ntin). The ratios 
for the A and K'^, shown by dotted-dashed and dashed lines, 
are taken from |1^ : Raa (b) and (c) are the ratios of central 
Au + Au (top 5%) to p -I- p and peripheral Au + Au (60- 
80%) to p + p, respectively. The values of Raa for charged 
hadrons are shown as open circles [25| . The width of the gray 
bands represent the uncertainties in the estimation of (Ntin) 
summed in quadrature with the normalization uncertainties 
of the spectra. Errors on the (j> data points are the statistical 
plus 15% systematic errors. 



ential cross-section in p-l-p is estimated as the NSD yield 
times (Jnsd, measured as 30.0 ± 3.5 mb, with a small 
correction, determined from Pythia calculations, of 1.05 
at pt=0.4 GeV/c and unity above 1.2 GeV/c J25|. Rcp 
is the ratio of the yields between two Au-|-Au centrality 
classes, scaled by {Nbin}- The Rcp (FigureO^a)) for the 
meson at moderate pt (1.5 < pt < 4 GeV/c) is sup- 
pressed relative to the binary collision scaling (dashed 
horizontal line at unity). 

A comparison of the Rcp for the </>, Kg and A is shown 
in figure O (a). Both statistical and systematic errors are 
included in the figure. The ratio Raa for central (top 
5%) and peripheral (60-80%) Au + Au data are shown 
in Figure |31 (b) and (c) , respectively. Raa for charged 
hadrons 25] is also shown as a reference. The charged 
hadron and 4> peripheral Raa both go above the binary 
scaling limit, but are consistent with unity within the 
systematic uncertainties. The (f) central Raa approaches 
unity and point to point is higher than Rcp- With the 
systematic uncertainty on the normalization of the ratio, 
however, both Raa and Rcp are consistent. Note that 
a Raa ratio that is higher than the Rcp ratio would be 
consistent with OZI suppression of </> production in p + p 
m, 1^ 1^ and/or strangeness enhancement in Au + Au 



collisions. A measurement of Raa vs. system size may 
be sensitive to the system size at which OZI becomes 
irrelevant to (f) production. 

The (f) Rcp result is consistent with a partonic recom- 
bination scenario 30, 32, 33]. In these models, the cen- 
trality dependence of the yield at intermediate pt depends 
more strongly on the number of constituent quarks than 
on the particle mass. Further higher statistical data for 
the (f) are needed to draw a conclusion. 

In summary, STAR has measured (p meson production 
in yfs^ = 200 GeV Au + Au and NSD p + p collisions 
at RHIC. The (f)/K^ yield ratios from e+e^, p + p and 
A + A collisions over a broad range of collision energy 
above the cj) production threshold are remarkably close to 
each other, cj) production, when scaled by the number of 
binary collisions, is suppressed with respect to peripheral 
collisions in central Au + Au collisions. The lack of a 
significant centrality dependence of the (/>/if ~ ratio and 
the values of {pt) effectively rule out kaon coalescence as 
a dominant production channel for the cj) at this energy. 
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